Electron Spin Resonance of defects in the Haldane System Y 2 BaNiOs 
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We calculate the electron paramagnetic resonance (EPR) spectra of the antiferromagnetic spin-1 
chain compound Y2BaNii_ a: Mg :E 05 for different values of x and temperature T much lower than the 
Haldane gap (~ 100K). The low-energy spectrum of an anisotropic Heisenberg Hamiltonian, with all 
parameters determined from experiment, has been solved using DMRG. The observed EPR spectra 
are quantitatively reproduced by this model. The presence of end-chain S = 1/2 states is clearly 
observed as the main peak in the spectrum and the remaining structure is completely understood. 



PACS numbers: 75.40 Cx, 75.10Jm, 75.40 Mg. 

A great deal of interest in one dimensional (ID) Heisen- 
berg chains with nearest-neighbor antiferromagnetic ex- 
change coupling, J, has been originated by Haldane's 
conjecture that integer- valued spin chains would exhibit 
a gap in the spin- wave excitation spectrum An in- 
teresting model of the underlying wave function for the 
Haldane state has been proposed by Affleck et al. 
They have shown that the exact ground state of the 
Hamiltonian ^iS -Si+i + ' Si+i) 2 /3 is a valence-bond- 
solid (VBS) state. This state in an open chain has two 
unpaired 5 = 1/2 spins, one at each end. In agree- 
ment with this, exact diagonalization || of finite open 
Heisenberg chains (without the biquadratic term) shows 
that the four lowest-lying states are a triplet and singlet 
whose energy separation approaches zero exponentially 
with increasing length. Monte Carlo B and density- 
matrix renormalization-group (DMRG) la] studies clearly 
show the presence of 5 = 1/2 end states. 

This picture was supported by EPR measurements 
of [Ni(C 2 H 8 N 2 )2(N02)]C104 (NENP) doped with non- 
magnetic ions j7j , where resonances corresponding to the 
fractional spin 5 = 1/2 states at the "open" ends of the 
5 = 1 Ni chains were observed. Similar measurements 
for doping with magnetic ions are also consistent with 
5 = 1/2 end states [Q. This effect, if robust, would cor- 
respond to the only instance in magnetism where a low 
energy collective excitation has no classical analog. 

However, Ramirez et al. || also tested the presence 
of free 5=1/2 states by studying the specific heat of 
non-magnetic defects in Y2BaNi05, with magnetic fields 
up to 9T and temperatures down to 0.2K. They found 
that the shape and magnitude of the Schottky anomaly 
associated with the defects in Y2BaNii_ a; Zn. I ;05 are bet- 
ter described by a simple model involving spin-1 exci- 
tations, instead of the 5 = 1/2 excitations of the VBS. 
In order to eliminate the apparent discrepancy with the 
EPR measurements in NENP, Ramirez et al. ||] pointed 
out the possibility that a small fraction of ethylene di- 
amine complexes in NENP acquire charge at structural 
defects induced by Zn doping. As it is not uncommon 



that structural defects induce a paramagnetic behavior 
in organic compounds, they propose that similar EPR 
measurements should be performed on Y2BaNii_ a; A a ;05 
(A a nonmagmetic ion) where x is more easily calibrated 
M. EPR is an appropriate technique to determine the 
existence of S=l/2 end states because it allows to distin- 
guish between S=l/2 and S=l spins in the presence of 
spatial anisotropy. 

In a recent paper || , a precise fit of the above men- 
tioned specific heat measurements was done, solving the 
low-energy spectrum of an anisotropic Heisenberg Hamil- 
tonian for Y2BaNii_ 2; Zn 2 ;05. It was shown that there 
is no contradiction between EPR in NENP |?]] and spe- 
cific heat measurements in Y 2 BaNii_. E Zn ;r C)5 ||. The 
results supported the existence of 5 = 1/2 excitations 
for sufficiently long chains and clearly indicated that the 
anisotropy plays a very important role in the low tem- 
perature properties. However, EPR experiments would 
unambiguously detect the presence of such excitations 
and definitely confirm the validity of the theory. 

Recently, Saylor et al |b|] have measured the EPR 
spectra of Y2BaNii_ :c Mg : i;05 at temperature T — 2K for 
different values of Mg concentration. The spectra show a 
prominent main peak, which can be associated with free 
S=l/2 spins, surrounded by some secondary peaks. The 
appearance of the latter is not related with the interac- 
tion between the S=l/2 end states belonging to different 
chains. Thus, it is necessary to give an explanation for 
the origin of these other resonances. 

In this paper, we calculate the EPR spectra of 
Y 2 BaNii_ a; Mg a; C)5 at T=2K for different values of x 
(x = 0.002, 0.006, and 0.02) and compare them with 
experiment. Mg and Zn are non-magnetic impurities, 
so the substitution of these ions for Ni provides a sim- 
ple break in the chain. Due to the high sensitivity of 
EPR experiments, it is optimal for discerning between 
long chains (where the 5=1/2 end spins are nearly in- 
dependent) and the shorter chains (where the 5=1/2 
objects are overlapping). We use the same Hamilto- 
nian as in Ref. p], with the same parameters as those 
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taken there from inelastic neutron scattering [ |l 0|Jl l| . The 
calculated EPR spectra show several contributions from 
chains with different lengths which are quantitatively re- 
produced. These contributions give rise to the secondary 
peaks mentioned above. The location of these peaks are 
related to the effect of anisotropy on chains with differ- 
ent lengths. In this way, these secondary peaks constitute 
clear evidence that the observed spectrum is originated 
by the end states of the ID MO5 chains and not by para- 
magnetic impurities. 

Y2BaNi05 has an orthorhombic crystal structure with 
the Ni 2+ (S =1) ions arranged in linear chains with a 
nearest-neighbor antiferromagnetic superexchange cou- 
pling J. The interchain coupling J±_ is at least three 
orders of magnitude weaker, making this compound an 
ideal one-dimensional antiferromagnetic chain. While 
each Ni atom is surrounded by six O atoms in near oc- 
tahedron coordination, the true site symmetry is D^h- 
Neglecting J± (the effect of which will be discussed 
later), the appropriate Hamiltonian for a chain segment 
of length N between two Mg impurities is [|l0| — |l2|| : 

H = J2i JSi ■ St+i + D{Slf + E[(Sf) 2 - (Sf) 2 }} 

i 

va 

where z is along the chain axis, St is the total spin and 
g av is the gyromagnetic tensor . Recent estimates based 
on fits of the Haldane gaps (in x, y, and z directions) 
measured by inelastic neutron scattering , indicate J ~ 

280K , D .038J , and E .0127 J jlfljilft. As 

Y 2 BaNi0 5 has a Haldane gap - 100K @Jll), the spin 
wave contribution to the EPR is negligible at 2K. In 
this temperature range, EPR is dominated by the effect 
of the defects. For this reason, it is necessary to solve 
accurately the low energy spectra (ui < 10K) of H for 
all values of ./V to be able to explain the observed EPR 
experiment. 

The lowest energy states of Eq.(Q) are a spin singlet, 
which we denote as |0), and three components of a split 
spin triplet (|15 t z )). The energy of |11) is the same as 
that of 1 1 — 1 ) due to the time- reversal symmetry of H. 
By means of the DMRG method, we have calculated the 
energy of these four states for all N < 40 and E = 0. 
The difference of energy e between any pair of eigen- 
states decays exponentially to zero with increasing N. 
This behavior allows us to extrapolate the energies to all 
values of N > 40, and demonstrate that the S —1/2 spins 
at the end of the chains are also asymptotically free in 
the presence of anisotropy. This issue is easy to under- 
stand considering that S =1/2 spins cannot be affected 
by anisotropy due to time-reversal symmetry. At this 
level, it is important to remark that the two end S =1/2 
spin excitations, which are localized at the end sites in 
the VBS, have a finite localization length I ~ 6 sites for 



the pure Heisenberg model |3| (5|. Therefore, while they 
are nearly free for large (N >> I) open chains, the inter- 
action between them is considerable when the length of 
the chain N is comparable to 21. This interaction splits 
the two S =1/2 states into a singlet (ground state for 
even N) and a triplet (ground state for odd N) ]I3] . 

The difference between any two energies of the above 
mentioned low-energy states is linear in D, and the 
quadratic corrections are negligible [jl2). This widely jus- 
tifies the validity of perturbation theory to first order in 
D. Then, by symmetry we can also include the term 

J2i E i( S i) 2 - ( s ?) 2 } t0 first OTder - Thus we find the 
following low-energy effective Hamiltonian including the 

triplet \1S Z ) and the singlet state |0): 



H eff = E (N) + (Ja(N) + Df3{N)) |0) <0| 
+E 1 (N)(Sl-Sl)-n B Y J B 
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Dl{N)Sl 
(2) 



where Eq(N), a(N), f3(N) and 'y(N) are functions of the 
chain length N, determined from the DMRG data fitting 
exactly the four lowest energy levels for E = 0. The 
numbers a(N), f3(N) and 'y(N) are represented in Fig. 
1. The validity of the last term of Eq.(||) has been verified 
explicitly by calculating the matrix elements of S^~, and 
Sf for all chains. H e ff determines the thermodynamics 
of the system at temperatures well below the Haldane 
gap. 
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Figure 1 

FIG. 1. Numerical coefficients of the effective Hamiltonian 
Eq. (Q) as functions of the chain length. 

Assuming a random distribution of defects, the prob- 
ability to find a Mg atom at any site i followed by N Ni 
sites and another Mg atom at site i + N + 1 is clearly 
x 2 (l — x) , where x is the concentration of Mg (missing 
S = 1 Ni spins). Thus, the probability per Mg impurity 
of finding a segment of N Ni atoms is x(l — x) . Then, 
the EPR signal intensity I(x, u, T) per impurity is: 
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I(x,u,T)=J2 x(l - x) N I N {u,T) 



(3) 



N=l 



where In(uj, T) is the EPR signal intensity of a segment 
of length TV described by H e f t . This intensity is given 
by the following expression: 



-0E a 



5{u-E u +E a ) 



(4) 



where \a) are the eigenstates of H (or H e jf) for the N- 
site chain, n is the direction (orthogonal to B) of the 
applied microwave magnetic field and (3 the inverse tem- 
perature. The delta distributions have been replaced by 
Lorentzian functions with a finite width to simulate the 
effect of interactions not included in Eq. (j^) and the ex- 
perimental resolution. 
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Figure 2 

Comparison between the experimental and the- 
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FIG. 2. 

oretical EPR spectra for three different concentrations of 
Mg, in the absence of interchain coupling. The values of 
g aa = 2.25A 9 bb = 2.177, and g cc = 2.165 are taken 
from Ref. [n3|. The parameters of the Hamiltonian (Eq . 
[l]) are taken from neutron scattering experiments [ ftC)| , [tl| : 
J = 280K, D = -0.039J, E = -0.127J. The arrows in- 
dicate the positions of the secondary peaks. 

In Fig. 2 we show the experimental data Jl^ | 
for I(x, T) with three different values of x (x = 
0.002, 0.006, and 0.02) at T=2K and compare them 
with the calculations. The uniform magnetic field B 
points in the x (or c) direction, while the direction of 
the microwave field is n=z (or a). The width and inten- 
sity of the lorentzians corresponding to each transition 
have been adjusted considering the main peak. These 
are the only adjustable parameters. The integrated in- 
tensity is determined by x. The most prominent feature 
of Fig. 2 is the main peak at around 8000G which is re- 
lated to very long chains for which the splitting induced 
by the anisotropy terms is much smaller than g cc [iBB. 



Surrounding the main peak there are some secondary 
peaks which correspond to the chains for which the mag- 
nitude of the anisotropy splitting is of the same order 
as g cc fiBB. In order to clarify the last statement, let 
us consider the simplest case of E = 0. The splitting 
A (AT) (between the S z = and the S z = ±1 states) in- 
duced by the anisotropy term DS^ (D < 0) for odd and 
even chains is represented schematically in Fig. 3. With 
respect to these latter states, a uniform magnetic field 
B along the x direction (perpendicular to the chains), 
does not shift the S x = Q state ((|11) - |1 - l))/y/2) but 
mixes 1 10} and -^(|11) + |1 — 1) giving rise to a split- 
ting between the lowest triplet states Ti and S x = 0: 
S(N) = A(iV)/2 + v /(A(AT)/2) 2 + {g cc /i B B) 2 . The appli- 
cation of a microwave magnetic field along the direction 
of the chains will induce transitions between these two 
states giving rise to resonant peaks when hu> = 6(N), 
being u> the microwave frequency. Then, chains larger 
than 40 sites (A(iV) << g cc fiBB) contribute to the main 
peak centered at htu = g cc ^sB. For chain lengths be- 
tween 15 and 30, A (AT) is of the same order as g cc ^sB 
giving rise to some secondary peaks at higher and lower 
energies than g cc \XBB. A(N) is negative for odd chains 
and positive for even ones, so the peaks on the left side 
with respect to the main peak come from even chains 
while the ones on the right side are produced by odd 
chains. In Fig. 2 c) we have also plotted the absorption 
curve corresponding to no anisotropy (D = E = 0). As 
is expected from the analysis given above, no feature is 
observed apart from the main peak. We also note that 
the weight of the secondary peaks is bigger for larger 
doping concentration (x — 0.02). This is due to the fact 
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FIG. 3. Diagram which shows the effect on the triplet 
states of an applied magnetic field B in the x direction in 
the presence of the anisotropy term Sf. 
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The calculated spectra are very sensitive to the values 
of D and E. However the best agreement is obtained 
using the experimental values given above @[n). In 
order to show the sensitiveness of the spectra to these 
parameters, in Fig. 4 we plot the spectrum for x = 0.02 
and two different sets of parameters. 
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Figure 4 

FIG. 4. same as Fig.2 for x = 0.02 and a) D = 0.039J, 
E = 0.015J; b) D = 0.030J, E = 0.0127J. 

The position of the different peaks is well reproduced 
within an error of 1000G^ IK (the spike at around 
8500G is a calibration feature). This error is of the same 
order of magnitude as the interactions neglected in our 
theory |l6| . It is important to note that there are no fit- 
ting parameters changing the positions of the secondary 
peaks. This agreement shows clearly that the secondary 
peaks of the EPR spectra are originated by the effect 
of spatial anisotropy on chains with lengths between 15 
and 40 sites. Due to experimental difficulties at very low 
fields, there is an artificial cutoff below 2500G. 

In Ref. B, we have included the interactions between 
chain segments in order to reproduce the specific heat 
measurements of Ref. || at zero magnetic field. Because 
of this it is important to determine how sensitive the cal- 
culated EPR spectra are to the inclusion of these interac- 
tions. The origin and the magnitude of this J' inter-chain 
interaction is explained in Refs. ||[l5]]. Since a calcula- 
tion of all possible J's and their effect on a 2D (or 3D) 
topology is a formidable task, we model the effect of J' 
considering a collection of chain segments with effective 
ID topology with exchange interaction J', coming from 
a uniform distribution between — J^iax an d ^max- We ob- 
tain that the inclusion of realistic values of J' has no 
effect on the position of the peaks. The only effect of 
J', as might be expected, is a broadening of the sec- 
ondary peaks and the corresponding small reduction of 
their maxima. As a consequence, the agreement between 
theory and experiment slightly improves for the peaks 
near 6000G and 12500G. 



In conclusion, by solving the low-energy spectrum of 
a Heisenberg Hamiltonian H which includes experimen- 
tal axial and planar anisotropy, we have reproduced 
in detail the low-temperature EPR data measured in 
Y2BaNii_ K Mn :E 05 for different concentrations of Mg. It 
is important to remark that there are no fitting parame- 
ters to reproduce the positions of the different peaks in the 
EPR spectra. These results confirm that the S =1/2 end 
chain excitations, which are asymptotically free for large 
chain segments, are experimentally observed. The fact 
that we have been able to fit all peaks with this model 
indicates that the 5=1/2 feature observed by EPR ex- 
periments is unambiguously an end-chain excitation, rul- 
ing out any other possibility. Summarizing, the EPR 
experiments can be interpreted in terms of end states of 
an anisotropic Heisenberg model, and the discrepancies 
about the existence of spin 1/2 excitations without clas- 
sical analog are removed. 
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